Noise in an insect outbreak model 
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We study the steady state properties of an insect (spruce budworm) 
outbreak model in the presence of Gaussian white noise. Based on the cor- 
responding Fokker-Planck equation the steady state solution of the proba- 
bility distribution function and its extrema have been investigated. It was 
found that fluctuations of the insect birth rate reduces the population of 
the insects while fluctuations of predation rate and the noise correlation can 
prevent the population of the insects from going into extinction. Noise in 
the model can induce a phase transition. 
Pacs numbers: 87. 10. +e, 05. 40. -a, 02. 50. Ey. 

I. INTRODUCTION 

Recently, Nonlinear dynamical systems with random noise have been paid much atten- 
tion both theoretically and experimentally. Phenomena such as noise-induced transitions, 
stochastic resonance, resonant activation, noise-induced spatial patterns are a few ex- 
amples of extensive investigations [1] [2] [3]. In most of these areas the noise affects 
the dynamics through system variables, i.e., the noise is multiplicative in nature [4,5]. 
The focal theme of these types investigation is the steady state properties of systems 
in which the fluctuations, generally applied from outside, are considered, independent of 
the system's characteristic dissipation. On the level of a Langevin-type description of 
a dynamical system, the presence of noise can change the dynamics of the system [6,7]. 
Noise processes have found applications in a broad range of studies, such as the steady 
state properties of a single mode laser [8], bistable kinetics [9], directed motion in spa- 
tially symmetric periodic potentials [10], stochastic resonance in linear systems [11], and 
steady state entropy production [12]. In this paper we study an insect outbreak model 
in the presence of the correlated noise, and show how noise can dynamically affect the 
ecosystem. 



II. AN INSECT OUTBREAK MODEL: SPRUCE BUDWORM 



A practical model which exhibits two positive linearly steady state populations is that 
for the spruce budworm, which can, with ferocous efficiency, defoliate the balsam fir. 
They are major problem in Canada [13], Ludwing [14] (1978) considered the budworm 
population dynamics to be governed by the equation 

f =r B N(l-^-)-P(N). (1) 

Here r# is the linear birth rate of the budworm and Kb is carrying capacity which is 
related to the density of the foliage available on trees. The P(N) term represents pre- 
dation, generated by birds. We take the form for P(N) suggested by Ludwing (1978), 
namely BN 2 /(A 2 + iV 2 ), where A is a positive constant and B represents the predation 
rate of the birds. The dynamics of N(t) is then governed by 



dN Ar/1 N . BN 2 

Before analysing the model we will express it in nondimensional terms. Here we in- 
troduce nondimensional quantities by 



x = = Ar B ,q = ^-,t = ^,(3 = B. (3) 



Upon on substitution (2) becomes 



^ = rx (i _ *) _ J?£_ ( 4 ) 
dr V q J 1 + x 2 



Now, if some environmental external disturbances make both the birth rate and the 
predation rate fluctuate, it is likely to affect r and (3 in the form of multiplicative noises 
that are connected through a correlation parameter A. As a result we have 

Where T(t) and £(t)are Gaussian white noises with the following properties: 

<r(f)> = m) = o, (6) 

(T(tns)) = 2D6(t-s), (7) 
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{Z(t)Z(s))=2a6{t-s), 



(8) 



mm) = (mm) = 2Av / ^^ - *). (9) 

Where D and a are the strength of noise T(t) and respectively, and A denotes 
the degree of correlation between T(t) and with < A < 1. 

We can derive the corresponding Fokker-planck equation for evolution of Steady Prob- 
ability Distribution (SPDF) based on Eq. (5)-Eq.(9). The equation is as follows [15] 

dP(x,t) dA(x)P(x,t) d 2 B(x)P(x,t) 

di ~ dx + dx~ 2 ' ( ' 

where 

A(x) = h(x) + Dg^g'^x) + \\j r Dag 1 (x)g' 2 (x) + \VDag' 1 (x)g 2 (x) + ag 2 (x)g' 2 (x), (11) 

B(x) = Dgf(x) + 2\VD^ gi (x)g 2 (x) + ag 2 2 (x). (12) 

Here h(x) = rx{\ - | ) - g x (x) = x(l - x/q), g 2 {x) = The steady 

probability distribution of the Fokker-Planck equation is given by [15] 

/ % N r f x A(x') , 

PJx) = W) eM J W) dx] - (13) 

Where Nq is the normalization constant. 



III. THE FLUCTUATION OF BIRTH RATE IN THE MODEL 

If we only consider the fluctuation of birth rate on the model, namely a = and A = 0, 
we can get 

- )~^-^ + Dx 1-- 1- — 



= rx(l - -) - + --)(!-—), (14) 



£(a;) = L>b(l - -)] 2 . (15) 

From Eq. (13), using the forms of A(x) and B(x), we get the following integral forms of 
the SPDF [16]. 

P^x) = w ^e X p[—]. (16) 

Here 
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X , 



9i(x) = x(l - -), 



(17) 



f(x) = r In 



x 



+ 



Pq 2 



x — q 1 + q 



(arctanx H ) + 



2[3q 3 



(In 



x — q 



x-q' (l + g 2 ) 2V 'v/TTx 



q arctanx) 



(18) 



The extrema of the SPDF are calculated using the condition A(x) — B ' (x) = 



r(l--)--f^- J D(l--)(l--) = 0. 
q 1 + x l q q 



(19) 



As for the parameters of the above equations we adopt r = 1.0, q = 10.0, (3 = 2.0, The 
results are represented in Fig.l-Fig.2 
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FIG. 1. Plot of the extrema of SPDF as a function of x for different noise strength values: 
D = 0.0, 0.1, 0.3, 0.7,using (3 = 2.0, and q = 10.0. 

In Fig.l, for zero noise strength, the curve at r = 1.0 gives only one x value, but as 
the value of D increases, the curve changes. The curve shows three x values at r = 1.0. 
From the figure we can see that the noise strength can change the state of the system. 



In Fig. 2 we show the effect of the noise strength D on the SPDF. For a small value 
of D, the SPDF shows a single peak region, which changes for larger values of D (see 
Fig. 2). As the value of D increases, the peak for the larger values of x decreases. At 
the same time, for small values of x, a new peak appears. It can be said that the noise 
causes the system become from two states to one state, namely noise can induce a phase 
transition. On the other hand, since x, D give the relative budworm population and the 
fluctuation of the budworm's birth rate, respectively, we can see from the figure the effect 
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the fluctuation of the birth rate makes on the growth of budworms, they can even make 
the budworms go into extinction. 




FIG. 2. Plot of P s t(x) (denotes the probability) against x for different noise strength values: 
D = 0.04, 0.1, 0.3, 0.6, using r = 1.0, (3 = 2.0, and q = 10.0. 

IV. THE FLUCTUATIONS OF THE PREDATION RATE IN THE MODEL 

On the other hand, if only the fluctuation of predation rate is investigated, namely 
D = and A = we can get the stationary probability distribution function similar to 
Eq.(14) 



Pst{x) 



\92{X) 



exp[ 



/(*)• 



(20) 



Where 



92{x) 



1 + x 2 ' 



(21) 



■x 3 + ^x 2 -(/3+ —)x + 2r\nx + (,'> 



)x x 

2 



(22) 



3q 2 q q 

In order to discuss the effect of the fluctuation of the predation rate on the system we 
adopt r = 1.0, q = 10.0,(3 = 2.26, the results are shown in Fig.3. 



Fig. 3 shows the effect of the noise strength a on the SPDF. For a small value of a, 
the SPDF shows the typical bistable region (see Fig.3) which vanishes for large values of 
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a. As the value of a increases the peak at the small x value decreases drastically, while 
for a large a value the curve has a single peak at the large x value . Since x denotes the 
budworm relative population, it is clear from Fig. 3 that, with an increase in the fluctua- 
tions of the predation rate, the budworm recover from going into extinction. 
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FIG. 3. Plot of P s t(x) (denotes the probability) against x for different noise strength values: 
a = 0.04, 0.10, 0.70 using r = 1.0, (5 = 2.26, and q = 10.0. 



V. EFFECT OF NOISE CORRELATION IN THE SYSTEM 

Since the two type of fluctuations have the common origin, (environmental external 
disturbance), we will consider the effect of the correlation between the birth rate fluc- 
tuation and the predation rate fluctuation. Based on Eq. (10)-Eq. (13) the stationary 
probabilities distribution function is represented by Fig. 4 with r = 1.0, (3 = 2.0, D = 0.3, 
a = 0.3. 

Fig. 4 gives the effect of noise correlation parameter A on SPDF. For a small A value 
the SPDF shows a typical bistable region which will vanish for the large value of A. As 
the correlation parameter increases the height of the peak on the small x value decreases, 
while the height of the peak on the large x value increases, namely probabilities flow from 
the small x value to the large x value. Since x denotes the budworm relative population, 
it is clear from Fig. 4 that the correlation of the noises is a advantageous factor for the 
growth of the budworm. 
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FIG. 4. Plot of P s t{x) (denotes the probability) against x for different correlation parameter 
values: A = 0, 0.4, 0.70, 1.0 using r = 1.0, (3 = 2.0 ,q = 10.0, D = 0.3, and a = 0.3. 

VI. SUMMARY 

In this paper, we studied the steady state properties of an insect outbreak model in 
the presence of the correlated noise. Birth rate fluctuation, predation rate fluctuation 
and the correlation of the noises are investigated. Birth rate fluctuation work against the 
growth of the budworm population, while both fluctuation of the predation rate and the 
correlation of the noises can prevent the budworm population from going into extinction. 
On the other hand, noise can dynamically affect the ecosystem. The noise from the birth 
rate fluctuation can make the system change from a single steady state to a bistable state, 
while noise from the predation rate fluctuation induces the system change from a bistable 
state to a single steady state. Thus in the insect outbreak model, noise from the fluctua- 
tion of the parameters can induce a phase transition. This viewpoint is completely novel 
in the traditional viewpoint it is the stochastic force that disturbs the phase transition. 
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